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The relativistic correction of the AdS/CFT implied heavy quark potential is examined within 
the framework of the potential model. For the typical range of the coupling strength appropri- 
ate to heavy-ion collisions, we find the correction is significant in size and lowers the dissociation 
temperature of quarkonia. 

! PACS numbers: 

o 

O ■ 1. INTRODUCTION 

O: 

jy-j ' The phase structure of quantum chromo dynamics (QCD) remains an active field of researches. At sufficiently 
fSl , high temperature, hadronic matter will evolve into a quark-gluon plasma which has been explored experimentally by 
relativistic heavy ion collisions (RHIC). The quarkonium dissociation is an important signal of this transition [1]. 

The quarkonium are bound states composed of a heavy quark q and its antiparticle q. It is found that, the ground 
states and the excitation levels of quarkonium are very much smaller than the normal hadrons, and that they are very 
tightly bounded 01 • Theoretically, there are two approaches to study the quarkonium: Lattice QCD and potential 
^4 models iSj. From Lattice QCD, we can calculate the spectral function numerically via the quarkonium correlators and 
identify the quarkonium states with the resonance peaks|^-t8,]. The potential model relies on the small velocity(t; << 1) 
of the constituent quarks. By solving a non-relativistic Schroedinger equation with a temperature dependent effective 
potential, we can determine the energy levels and thereby the threshold temperature when the bound state dissolves 
|9l-[l4j. The potential model will be applied in the present work. 

AdS/CFT correspondence is a powerful tool to explore the strongly coupled A/" = 4 super Yang-Mills plasma. 
The equation of state and viscosity ratio, etc. extracted from the AdS/CFT show remarkable agreement with the 
. Lattice QCD or experimental data from the quark-gluon plasma (QGP) created via Relativistic Heavy Ion Collisions 
\^ ' (RHIC). It would be interesting to extend the comparison to a wide range of other quatities, for instance heavy quark 
. dissociation, which are calculable in both ways to assess whether the super Yang-Mills serves an important reference 
■ model of the QGP phase of QCD. This is the primary motivation of this paper. 
fS| '. In a- previous workTs], we examined the heavy quarkonium dissociation within the potential model with the 
T— I ■ AdS/CFT implied potential function (holographic potential). We found that the holographic potential can be ap- 
^ \ proximated by a truncated Coulomb potential to a great accuracy. With the typical values of the 't Hooft coupling 
• w^ . constant, A = \/ NcQym considered in the literature [l7| . 

5.5<A<67r, (1.1) 
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our dissociation temperatures are systematically lower, though not far from the lattice prediction. [29[ On the other 
hand, an estimate of the velocity of the constituent quarks inside the bound state indicates that the non-relativistic 
approximation may be marginal, especially for J/ 5*. This motivates us to examine the relativistic corrections of the 
holographic potential with the aid of a two-body Dirac equation(TBDE). 

While the holographic potential alone is sufficient in the nonrelativistic limit, it does not provide all information 
necessary for the relativistic corrections even to the order v'^ term. Except for the correction brought about by the 
relativistic kinetic energy, the spin-orbital coupling and the Darwin term depend on how the holographic potential 
is introduced in the two-body Dirac equation. In addition, gravity dual of spin dependent forces is not available in 
the literature. Therefore our result remains incomplete at this stage. We would like to comment that the same issues 
exists for the relativistic corrections of the heavy quark potential extracted from the lattice QCD simulations. 
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In the next section, the work reported in [15| will be reviewed and the dissociation temperature beyond the 
truncated Coulomb approximation is presented. The corrections to the dissocation temperature are computed in 
the section 3 through a Foldy-Wouthuysen (F.W.) transformation of a two-body Dirac Hamiltonian. The kinetic 
energy contribution and the contribution from the Darwin and the spin-orbit coupling are calculated seperately with 
the latter obtained simply by replacing the perturbative Coulomb potential in the two-body Dirac Hamiltonian with 
the holographic potential. The section 4 conclude the paper. 



2. THE HOLOGRAPHIC POTENTIAL MODEL 



In the conventional potential model of QCD, the non-relativistic wave function of a heavy quarkonium satisfies the 
Schroedinger equation 

h^V^ + t/(r, T)]^ = -Em, (2.1) 

where E{T) is the binding energy and U{r,T) is identified with the internal energy of a pair a static q and q in QGP 
and is related to the free energy F{r, T) via 

U{r,T) =-r'[±{l^)^^ (2.2) 

The free energy F{r,T) can be extracted from the expectation of a pair of Wilson loops operator according to: 



e 



-iF(r,T) _ 



tr < W^L+)W{L^) > 
tr < WHL+) >< WiL_) > 



(2.3) 



where L± stands for the Wilson loop running in Euclidean time direction at spatial coordinates (0,0, ±^r) and is 
closed with the periodicity /3 = i and 

M^(i±)=Pe"'^-±'"'^^^"\ (2.4) 

The spatial coordinates of L± are (0, 0, ±|r). The lattice QCD simulation of the expectation value (|2.3[) can be found 
if Ref.[l3,[ll. 

In case of super Yang-Mills, the holographic principle places the Wilson lines L± on the boundary (y — oo) of the 
5D AdS-Schwarzschild metric p^ : 

ds^ = TT'TYifdt' + dx^)+ ,J,,J y^ (2.5) 
where / = 1 — dx^ = dxf + dx\ + dx\ with the ansatz xi — X2 = Q and X3, a function of y. 



The free energy F(r, T) of the corresponding super Yang-Mills at large Nc and large 't Hooft coupling is proportional 
to the minimum area of the worldsheet in the AdS bulk bounded by and L_, and is given parametrically by ^19i.i20i]: 

F{r,T) = Tmin(/,0) 

_ 2q dy (2.6) 

where 



I = V\[r dyj ^4 - 1) + 1 - (2.7) 

and the parameter yc G (l,oo). Eliminating yc between (j2.6p and (j2.7p . we find that 

F{r,T)^~-^p)eipo-p) (2.8) 
where a = 0.2285\/A, p = nTr, pa = 0.7541 and $(p) is the screening factor. The corresponding internal energy is 

Uir,T) = -^[$(p) - p(|),. - )]^(Po - P) (2.9) 
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Td (A = 5.5) 


Td (A = Gtt) 




Exact 


truncated Coulomb 


Exact 


truncated Coulomb 


Is 


142 


143 


262 


265 


2s 


27 


27 


50 


50 


Ip 


31 


31 


57 


58 



TABLE I: It lists our melting temperature in MeV's for Is, 2s and Ip state using the exact potential and compared with the 
Coulomb potential case. We can find that this two results are very close to each other which confirmed that the Coulomb 
approximation is excellent. 



and will be substitute into the Schroedinger equation (|2.ip . 
The small p expansion of ^{p) is: 

Within the screening radius po, the first two terms of the series (|2.10p approximate the exact $ well as is shown in 
Fig.l. If we keep only the first two terms, the screening radius po — 0.7359 and U (r, T) becomes a truncated Coulomb 
potential 

u = --e{po^p) (2.11) 

r 

under the approximation. 

We define the dissociation temperature Td as the temperature when the binding energy falls to zero, i. e. E{Td) — 0, 
and the corresponding radial Schrodinger equation reads [15]: 



(PR 2dR Ml + l) 



dp^ p dp p^ 



V]R = Q (2.12) 



where the reduced potential, V = J^'^^ is dimensionless. 

The truncated Coulomb potential approximation was employed in fl5!| and the dissociation temperature of the 
bound state of Z-th partial wave and n-th radial quantum number is given by 

Td^'-^ (2.13) 

with Xni the n-th nonzero root (ascending order) of the Bessel function J2i{x). The corresponding radial wave function 
reads 

R{r) = -^J2i+i (xnuf^ ■ (2.14) 



VP 

for p < Po and R{r) = const. /r'+^ for p > po. 

In this work, we have calculated the dissociation temperature with the exact holographic potential (|2.9p . The 
comparison with that obtained from the truncated Coulomb potential in [l5| for J/^* is shown in Table I, where 
we choose m = 1.65 GeV for the mass of c quarks. From the comparison of these two results we confirmed that 
the truncated Coulomb approximation is a good approximation and we shall stay with the truncated Coulomb 
approximation for the rest of this paper. 



3. THE RELATIVISTIC CORRECTION OF THE HOLOGRAPHIC POTENTIAL 



As it is mentioned in the introduction, the velocity of the heavy quarks is not low enough so the relativistic correction 
may be significant, especially for J/vt. To explore this correction, one has to go beyond the Schroedinger equation 
(12. ip and switch to the two body Dirac eg nation [2l| - |2^ : 



(3.1) 
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In the center-of-mass frame, the Hamiltonian of the two body Dirac equation is: 

H = ai ■ p + /3i ■ m — a2 ■ P + 13-2 ■ m + U (3-2) 

where, di = a (8) a2 = ^ "8) a, /3i = /? i8) /, /32 = ^ <8) cJ, /? are usual 4x4 Dirac matrix, p — — iV and U is 
the interaction potential between the two particles. The Hamiltonian _ff is a 16 x 16 matrix. A quarkonium state 
corresponds to a bound state of H with the eigenvalue 2m — E{T), which goes to 2m at the dissociation temperature, 
i. e. E{T) = 0. Since we are interested in the leading order relativistic correction of the dissociation temperature for 
the quarkonium, we have to expand the Hamiltonian to the order w^. The sorting of the order in v follows from the 
rules that ^ ~ [/ ~ and V ^ ~ mv. Also the expectation values of cEi and a2 are of the order v. 
In analogous to the one body Foldy-Wouthuysen transformation we introduce the unitary operator 



where 



U = e'^'^e'^'^e'^^e'^' (3.3) 



(3.4) 
(3.5) 
(3.6) 
(3.7) 

with 

Oi 

The transformed Hamiltonian reads 
i/pw = UHU'^ 

ff^ 1 1 dU 

= il3i+l32)im+^-^) + U+-^V'U+-^ — {ai+a2)-L (3.8) 
2m 8m'' 4m^ 4TO^r dr 

where, higher order terms in v have been dropped. The details of the transformation are deferred to the appendix. 
The non-relativistic wave function, Vl* = \['sjS2(fi — r2) with subscript Si, S2 labelling the spin components of the two 
quarks, corresponding to the sector with /3i = /32 = 1 (This wave function can be expanded in the series of products 
of the orbital wave functions of the preceding section and the spin wave functions). We may stay within this sector 
for the 1st order perturbation of the v'^ term of (|3.8p with the effective Hamiltonian Hcs. = Hq + Hi, where 

-2 

Ha = 2m+ — + U 



Si = 




■Oi 


S2 = 


2m 


■ (-O2) 


S'l^ 




■o'l 


S'2^ 


2 m 


■ i-0'2) 


di -p 


S O2 = 


-- d2-p. 



m 



corresponds to the non-relativistic part, 
and 



1 4j^3 ^ 4^2 4jJj2j. ^j. ^ 1 + ^2) 

- ^ ' ^ ;V^U + -^'^{J^ - - S^), (3.9) 



4m'^ 4m^ Am^r dr 

is a relativistic correction. We have introduced the total spin S = + (^2) and the total angular momentum 

J = L + 5 in the last step. The contribution of Hi is somewhat like the that to the fine structure of the hydrogen 
atom, and we can see that the first term is the first order of the kinetic energy correction; the second term is Drawin 
term; and the third term is the spin orbit coupling which can be decomposed into the spin singlet and the spin triplet 
channels. 

Perturbatively, we may write E{T) — Eq{T) + 6E{T) and T = Tq + 5T, where Eq{T) is the non-relativistic binding 
energy in (12.11) . and Tg are the dissociation temperature given by (|2.13p . and 5E{T) and 5T are the corrections. 
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We have Eq{To) = 0. Expanding the dissociation condition E{T) = to the order w^, we obtain the formula for ST, 
i.e. 



where, 



5T = SiT + 62T =- 



6iE{To) + S2E{To) 

V QT ''-'0 



(3.10) 



and 



SiE{To) 
S2E{To) 



< -7-^ > 



< 



Am? 
1 

Am? 



V^U > + < 



dE^ 
dT 



1 dU 
Am^r dr 



T=To 



(3.11) 



(3.12) 



The average 



f d^ri;*(r)0(r)^/j(r) 

< O >= 



J d'^ril)*{r)^{r) 



(3.13) 



with '4!{r) the non-relativistic wave function. The reason for our separating the contribution from p^, SiT and that 
from the Darwin and spin-orbital terms, S2T is the uncertainty in the representation of the holographic potential in 
p.2p . which does not impact on the p"^ correction. We will come to this point in the next section. In the limit of zero 

binding energy, we find < >=< -^^W^ >■ For the truncated Coulomb potential, 



V^U = ATTaTT^T^S^{p)e{pQ - p) + 



-S'{p- po). 



(3.14) 



In terms of the radial wave function Ri{r) of '0(r'), 

„2 



< 



< 



< 



1 



-4 

Am^ 



> 



Am^ 

1 dU 
Am^r dr 



For the ns state, we find that 



dU 



SiT 



a 



V^U > = 



AmnT 
a 



Po 



Riipfdp 



> = 



Am^ 

a 
Am^ 

mi 



mm? - 2Ri{po)R'i{po)po - rKpo)} 

Po 



— R (p) + - — ^ 

p ^ ' Am^ 



R'ipo) 



Rf{Po)pl 



TTaT§ - 



Ampo 
SiT + S2T = 



Ampo 



For the np state, we can also get analytical expressions, which are more lengthy. 



(3.15) 



(3.16) 



The numerical values of the corrected temperature Tq + 5iT and Tq + 5iT + 52T in MeV's for Is, 2s and Ip states 
are listed in the Table II below. 



4. DISCUSSIONS 



In summary, we have explored the leading relativistic correction to the dissociation temperature of heavy quarko- 
nium state through a F.W.-like transformation of the two body Dirac Hamiltonian with the AdS/CFT implied 
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cc 


bb 




A = 5.5 


A = 67r 


A — 5.5 


A = 67r 


Is 


162.54 


387.54 


478.76 


1139.11 


2s 


29.15 


62.75 


85.67 


184.44 


Ip 


32.04 


62.14 


94.18 


182.66 






cc 


bb 




A = 5.5 


A = 67r 


A = 5.5 


A = 67r 


Isl 


130.79 


188.65 


385.63 


555.58 


Isl 


130.79 


188.65 


385.63 


555.58 


2sJ 


26.71 


48.16 


79.15 


142.59 


2s? 


26.71 


48.16 


79.15 


142.59 


Ipi 


31.53 


61.33 


93.54 


180.79 


IPo 


32.65 


68.48 


96.85 


201.80 




32.09 


64.90 


95.20 


191.30 




30.96 


57.76 


91.89 


170.29 



TABLE II: This lists the final results of the corrected temperature in MeV's, where the upper one corresponds to the results 
of To + SiT and the lower one corresponds to To + SiT + 52T. For the lower one, we wrote the states as nL^f^^ in the first 
column, where n is the main quantum number, L is the orbit angular momentum quantum number, S is the spin quantum 
number and J is the total angular momentum quantum number. Since there are no spin-orbit coupling term for ns states, 
the spin singlet and spin triplet are degenerate. However, for p state, adding the coupling term, there will have energy level 
splitting, so we can see the different results for different total angular momentum J. 

potential. Among the contributions we considered, the correction of the kinetic energy, being negative, enhances 
the binding but the Darwin term does the opposite and dominates. Consequently, the dissociation temperature of 
s-state is lowered, leaving the corrected values further below the lattice result. This disagreement can be attributed 
to the short screening length tq = about 0.25fm at T = 200 MeV, of the AdS/CFT potential and the sharp 
cutoff nature of the screening. In case of J/V', the magnitude of the correction ranges from 8% for A = 5.5 to 30% for 
A = Gtt, indicating significant relativistic towards the high end of the domain p.l3p . 

The potential model, though physically more transparent than spectral function approach, does not provide com- 
plete v"^ corrections with the holographic potential extracted from the Wilson loop alone. The same deficiency applies 
the relativistic correction based on the lattice heavy quark potential alone. As the Wilson loop for a nonAbelian 
theory involves multi-gluon excahnges, its form in the two-body Dirac Hamiltonian (j3.2p may not be adequate unless 
the single gluon exchange serves a reasonable approximation. Among the four sectors Pi — ±1 and /32 = ±1, only 
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Hi - -izr, + 732 V'C/+ + -r-T--j^{Si + <?2) • L + ^({^1 • V, {ai • V, [/_}} + {a, ■ V, {a, ■ V, [/_}} (4.2) 



two of them /3i = /32 = ±1 correspond to qq interacting via the holographic potential The other two sectors with /3i 
and P2 correspond to qq or qq and the interacting potential is unknown. A more general form of the interaction in 
p.2p without violating the charge conjugation symmetry is to replace U by 

(A++ + A__)J7 + (A+_ + K-+)U' = U++ pi(32U- (4.1) 

where the projection operator A±i_±i = ^^^^ ^^2^ ■, U' the potential between qq or qq and U± = . [/' = [/ in the 
previous section. Repeating the steps of FW transformation in the appendix, we find that the perturbing Hamiltonian 
p.9p is replaced by 

T-^ + 1—^^^+ + J-^ +$2)- L + — „ 

with {...} an anticommutator. This will modify the potential part of p.9p . Only the correction from the kinetic 
energy, term of (|3.9p is robust, which raise the dissociation temperature. 

In addition to the holographic potential considered in this work, there should be spin dependent ones that splits 
degeracies between spin singlets and spin triplets (e.g. between rjc and J/4')- The gravity dual of the latter are 
unknown in the literature. A first principle derivation of the spin-dependent forces associate them to the expectation 
value of Wilson loops with operator insertions, < tYW^v{L+)'^Wpx{L^) >, where 

= PF^,{x)e~'^L (4.3) 

with Ffj^^ the Yang-Mills field strength and x a point along L, W^i/(i) is obtained from the Eqs. (|2.4l) by small distortion 
of L at X. Within the AdS/CFT framework, it corresponds to the perturbation of the Nambu-Goto action of the 
world shee underlying the holographic potential under a small distortion of its boundary. It is a challenging boundary 
value problem and we hope to report our progress along this line in future. 

Finally, we would like to comment on a phenomenological formulation of the two-body Dirac equation 27] , which has 
been applied recently to the same problem addressed in this work [28] . It amounts to divide the heavy quark potential 
of Cornell type into a linearly confining term and a single gluon Coulomb term and to generate all spin-dependent 
forces by the latter. While it is legitimate in vacuum in the weak coupling because of the Lorentz invariance, a direct 
application to a medium beyond weak coupling remains to be justified, given different screening properties of the 
electric and magnetic gluons. 
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Appendix A 

In this appendix, we shall fill in the details of the F.W. transformation for the two body Dirac equation which we 
have done in section 4. 

Let us recall the Foldy-Wouthuysen transformation of the one body Dirac Hamiltonian, H = a ■ p -\- /3m -I- V . For 
a 4-component spinor with velocity w << 1 and positive energy, one can work in a representation where the upper 
two component corresponds to the non-relativistic limit, referred to as the large components, while the lower two 
components are suppressed by a power of v, referred to as the small components. An operator is even(odd) if it is 
diagonal(off-diagonal) with respect to large and small components. For example, /3 is even and a is odd. The Foldy- 
Wouthuysen transformation amounts to successive unitary transformations that push the odd operators to higher 
orders in v. 

The Foldy-Wouthuysen transformation can be easily generalized to the two body case, the Hilbert space of which 
is spanned by the direct products of two one body spinors. To the leading order relativistic correction, we stop at the 
terms ignoring all the higher order term. The Hamiltonian of the two body Dirac equation is: 

H = dip + j3im — a2P + + V (A.l) 

let: 

dtip = Oi 
a2P = O2 
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so, 

H =pim + Oi+ - O2 +V 

^ , ■' V , ' 

= Hi+H2 + V 

(A.2) 

where Hi corresponds to the first under-brace, and H2 corresponds to the second under-brace. And it is easy to proof 
the commutation relation that:[0i,02] = [Pi, (32] = [Oi,/32] = [02,/5i] = 0, and 01^02^ v. 
We do the first transformation: 
we select: 

51 =-^P.Oi 

52 =-i^H-0.) 

so, the transformed Hamiltonian turns into: H' = e^^^ e^^^He~^^^ e~^^^ . We calculate the mid over-brace first. 

^ V ' 

V V ^ 

(A.3) 

also, Hi corresponds to the first under-brace, H2 corresponds to the second, and V corresponds to the third. 

Ih =Hi + [iSi,Hi] + ^[iSi, [iSi,Hi]] + ^[iSi, [iSi, [iSi,Hi]]] + ^[iSi, [iSi, [iSi, [iSi,Hi]]]] + ■■■ (A.4) 

[iSi,Hi] = -Oi + -piOl 
m 

^[^Su[^Si,Hi]] = -±-PiOl-J-,Ol 

l^iiSi, [iSi, [iSi,Hi]]] = J-^Ol - ^^PiOi 

^[iSi, [iSi, [iSi, [iSi,Hi]]]] ~ ^I3i0t 

TTi =/3i- + i/3iO?-^/3iO^-^03 + 0(.^) (A.5) 

where v is the velocity of the particles. 

H'2 =H2 + [iSi,H2] + ^[iSi, [iSi, H2]] + ^[iSi, [iSi, [iSi, H2]]] + ^[iSi, [iSi, [iSi, [iSi, H2]]]] + ■■■ (A.6) 

[iSi,H2] = ^[tSi, [iSi,H2]] = ^[tSi, [iSi, [iSi,Hi]]] = • • • = 
H2 =H2 = 132m - O2 

so, 

= Am + ^2m - O2 + y + l^PiOl - ^[0„ [Oi,V]] - ^/3.0^ + ^A[0„ V] - ^O? + 0(.^) 
= H (A.8) 
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that we mark e^^^He as H for convenience. So, 



= H+ [iS2, H] + ^ [iS2, [iS2, H]] + 1 [iS2, [iS2, [iS2, H]]] + 1 [iS2, [iS2, [iS2, [iS2, H]]]] + ... (A.9) 

[iS2,H] =02 + ^P20l - l-p^[02,V] - ^/3l/32[02, [Ol,F]] 

^^^'^^ [^^2,^]] = -^P.Ol + - ^[02, [02, V]] 

^[iS2, [iS2, [iS2M\ = --^201 - ^^20| 

^[iS2, [lS2, [lS2, [iS2M]] - ^^20^ 

H' = Pim + /32m + V 

-^/3i/32[02,[Oi,y]] + 0(t;^) (A.IO) 



since, 



H' = + V + ±^.0l - ^[0„ [0„ y]] - ^/3iOf + ^ ^[0, F] - ^O?) ^ mar. : ^ 



1 . „o 1 _ _ 1 . 1 „ _ 1 

71 + ; 

1 



+(^2m + -/^.O,^ - —,[02, [02, V]] - —,P20t - -HO., V\ + ^ mar. : 



4^2-/3i/32[02,[Oi,y]] (A.ll) 

" V ' ^ V ' 

= /3im + y/ + O'l (A.12) 
H', = /32m + ±P20l - ^[02, [02,V]] - - (^^2[Q2, V] - ^Q^) 

= ^2m + + 0'2 (A. 13) 

H' =H[+H!,-^p,p2[02,[OuV]] (A.14) 

the first under-braces upside in the expression of -ff{ (iJa) correspond to V{{V2), and the second correspond to 0[{0'2) 
for convenience. 

Then, we do the second transformation: 
We select: 

and 0[ ~ ~ z;3 
so, 

= e''52 e»s^;^^e-'^i e"'"^^ +6*^= gisl^/g-isl ^-is'^ _ ^^is'^^is'^ f^^f^^p^^^ [Oi,V]]e-'^'^e-'^'^ 

= e'^'^Hie-'^'^ + e'^'^Hle-'^'^ - -^e'^^e''^'^PiP2[02, [Oi, F]]e-*^ie-'^^ (A.15) 
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where, the first under-brace upside correspond to H[, and the second correspond to H2, and we can directly apply 
the results of the first transformation 



= /3im + F/ + ^/3i[0;,F/] (A.16) 



= I3,m + V( + 1-(3,[0[,V;] - ±-l32[0'2,V{] (A.17) 
=H!, + [zS[,H^] + 0{v') 

= (3,m + Vi -0', + ^[0[,V^ - ^[O'M] (A.18) 



= + - ^p^p'M + ^fii[0'i,Vi] - ^p,[0'M (A.19) 



Then, we see the last term in H": e*^^ e^^'^ [d, V^]]e-'^i e"*^^ 

we do the calculation of the over-brace first. 

e'^^/3i/32[02, [Oi,V]]e~'^'^ = /3,/32[02, [Oi,V]] + [tS[, (3,/32[02, [Oi,V]]] + 0{v'') 

= (/3i/32[02, [Oi,F]] + ^[p,[0,,V], [O2, [0,,V]]]) =^ mark : A (A.20) 

gis^^g-»s^ =A + 0(w^) (A.21) 

^ f3if32[02AOuV]] (A.22) 



considering Oi, O2 ^ v, O'l, O2 



H" = /3im + V; + /32m + - ^/3i/32[02, [Oi,V]] 

= + /^^- + ^ + i/'iO? + ±P20l - [O, F]] - ^[02, [O2, VI 

1 . 1 . ^4 1 



8m3^^^i - ^/?20| - ^/3i/32[02, [Oi,F]] (A.23) 
smce, = Oi = f, Of ^ ^ 

[Oi,[Oi,l/]] = -VV--^s'iL 
r or 

[02,[02,V]]^-W^V--^f2L 

r or 

[02,[Oi,V]] = -aua2j\7,\/jV 

+ + E2)i + ^/3i/32aH«2,V.V,F (A.24) 

The last term in (jA.24[) . though of the order of v"^, is a direct product of two odd operators and therefore does not 
contribute to the first order perturbation considered in this paper. 
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